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We analyze a model for spin squeezing based on the so- 
called counter-twisting Hamiltonian, including the effects of 
dissipation and finite system size. We discuss the condi- 
tions under which the Heisenberg limit, i.e. phase sensitivity 
oc 1/N, can be achieved. A specific implementation of this 
model based on atom-atom interactions via quantized photon 
exchange is presented in detail. The resulting excitation cor- 
responds to the creation of spin-flipped atòmic pairs and can 
be used for fast generation of entangled atòmic ensembles, 
spin squeezing and applications in quantum information pro- 
cessing. The conditions for achieving strong spin squeezing 
with this mechanism are also analyzed. 

PACS numbers 03.67.-a, 42.50.-p, 42.50.Gy 



I. INTRODUCTION 

Interacting quantum systems that start in uncorrelated 
states generally evolve towards entangled states due to 
quantum correlations building up in time. These corre- 
lations and the form they take depend crucially on the 
interaction that gives rise to them. For example in para- 
metric down conversion or in the optical parametric os- 
cillator (OPO) pairs of photons can be created in distinct 
modes of the electromagnètic field. The fact that pairs of 
photons are generated leads to quantum correlations be- 
tween the two modes. Since each mode is described by a 
harmònic oscillator, one can think of the state of the field 
as the quantum state of two fictitious partides in har- 
mònic oscillator potentials. The quantum correlations 
correspond to e.g. the positions of the partides being 
strongly correlated, in the ideal case A(Ai — X2) 2 — > 
and their momenta being anticorrelated A(Pi+P2) 2 — > 0. 
For the electromagnètic field modes, the position and mo- 
menta correspond to quadratures of the field modes and 
it is between these that correlations are produced [jïl[|. 
These correlations are essential to quantum communi- 
cation e.g. quantum teleportation of information from 
one location to another || . Entanglcment is also crucial 
for many schemes in quantum cryptography and for long- 
distance quantum communication through lossy channels 

Since the mechanism for producing correlations in elec- 
tromagnètic field modes is at the fundamental level so 
simple (photons created in pairs) it is natural to wonder 



if such a simple mechanism may lead to entanglement 
of àtoms interacting in a similar manner. In complete 
analogy to the OPA mechanism, a process that transfers 
pairs of àtoms from their ground state to two well de- 
fined final states also gives rise to quantum correlations 
between àtoms. When a collection of N two level àtoms 
is thought of as an ensemble of effective spin ^ partides 
with total pseudo-angular momentum J = N/2, it turns 
out H that the quantum correlations produced by an 
interaction that transfers àtoms in pairs from the lower 
state to the upper state shows up as reduced fluctuations 
in a component of the angular momentum e.g. A J 2 — ► 0. 
We will discuss entanglement of àtoms with one another 
in an atòmic ensemble for which an effective interaction 
leads to the transfer of àtoms in pairs to well defined final 
states and we will use the concept of spin squeezing to 
quantify the amount of quantum correlations produced 
in such a case. As for squeezed states of light, deco- 
herence mechanisms and dissipation are acting in such 
a way as to destroy or limit the amount of squeezing 
achievable in practice. We also analyze the influence of 
such dissipation mechanisms and find relations between 
the spin squeezing interaction rate, the dissipation rate 
and the amount of squeezing achievable in the presence 
of damping mechanisms. The coherent control of the 
dynamical evolution of complex systems such as atòmic 
ensembles may lead to the production of entangled non- 
classical states such as spin squeezed states || (analogous 
to squeezed states of light Hj) an d correlated collective 
atòmic modes (similar to twin photons generated by a 
non-degenerate OPO). 

The main result of this paper is that for a collection 
of N àtoms with average single atom nonlinearity x (two 
atom interaction energy) and with single atom loss rate 
r, the condition for achieving some spin squeezing is that 
Nx > r. In order to achieve reduction of uncertainty in 
say J x compared to the uncertainty in the Bloch state 
|J = N/2,J Z = N/2) for which (AJ X ) 2 = N/A by an 
amount s (i.e. (AJ X ) 2 = N/(4s)) with 1 < s < N, one 
requires that Nx > sT and the interaction time needed 
seales as t ~ (log s) / {Nx) while the maximum number of 
àtoms than can be lost without destroying the squeezing 
seales as AiV ~ {N/s) logs. To achieve Heisenberg lim- 
ited precision (i.e. maximum spin squeezing s ~ N), one 
needs a large single atom nonlinearity x > T. This mcans 
that the interaction time needed to achieve this strongly 
correlated state is t ~ (log N)/(Nx) and the maximum 
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number of àtoms that can be lost without compromising 
this optimal squeezing is AN ~ logiV i.e. a very small 
number of àtoms lost may prevent reaching the Heisen- 
berg limit. This analysis remains vàlid and agrees with a 
spccific implemcntation based on an cffectivc atom-atom 
interaction via quantizcd photon exchange in a cavity, for 
which the decoherence mechanism corresponds to spon- 
tancous emission and leakage of photons from the cavity. 

The possibility of cohercntly controlling interacting 
quantum systems has lead to many new developments 
in the field of quantum information science ||. These 
are expected to have an impact in a broad area rang- 
ing from quantum computation and quantum communi- 
cation [Q to precision measurements || and controlled 
modeling of complex quantum phenomena |)|. Entan- 
gled systems realized in the laboratory range in size from 
few qbits [0, to macroscopic ensemble of partides [ïl·|. 
Controllable coherent interactions between àtoms rfL2 Jl3| 



may also open the way for modclling of complex quan- 
tum phenomena such as quantum phase transitions || in 
which quantum correlations play a crucial role. 

Entanglement of a single atòmic ensemble, i.e. quan- 
tum correlations between àtoms in the same ensemble, 
has been shown to be potentially very useful in the field of 
precision measurements Q . Certain types of interactions 
between àtoms lead to entanglement and spin squeezing, 
characterized by reduced variance in an observable and 
increased fluctuations in the canonically conjugate ob- 
servable. This rcduction of fluctuations directly trans- 
lates into an improved accuracy for measurements sensi- 
tive to that observable. A typical figure of mèrit for spin 
squeezed states is the phase accuracy 5<p on estimating 
accumulated dynamical phase in the Ramsey interfero- 
metric experiment. With all experimental uncertainties 
controlled below this noise level, the dominant source 
of noise in such experiments is the "quantum projection 
noise" j^] associated with e.g. the noise in measurements 
of the x-componcnt of the spin of an ensemble of two- 
level àtoms (effective spin i) all prepared in the lower 
level (the | J.) state). This noise leads to a lower limit on 
phase accuracy S(f> — 1/vN called the Standard quan- 
tum limit (SQL), where N is the number of àtoms in the 
ensemble. The Heisenberg uncertainty principle however 
allows for phase accuracies consistent with the bàsic prin- 
cipies of quantum mechanics that are as low as 5<j) = 1/N, 
called the Heisenberg limit. 

We also discuss in more detail a technique |Ï4| based 
on a resonantly enhanced nonlinear process involving Ra- 
man scattering into a "slow" optical mode [^5[ , which cre- 
ates a pair of spin-flipped atom and slowly propagating 
coupled excitation of light and matter (dark-state polari- 
ton) . When the group velocity of the polariton is reduced 
to zero |nï,0|, this results in pairs of spin flipped àtoms. 
The dark-state polariton can be easily converted into cor- 
responding states of photon wavepackets "on demand" 
[Q7f , which makes the present approach most suitable for 



implementing protocols in quantum information process- 
ing that require a combination of deterministic sources of 
entangled states and long-lived quantum memory [f|jïl| . 

This paper is divided into V sections. In Section II, 
we discuss Ramsey spectroscopy and the use of spin- 
squeezed states in precision measurements. In particular 
we analyze the situation where N two level àtoms with 
levels \g) and |e) are prepared in a correlated state and 
subsequently probed by separated fields of frequency ui in 
the Ramsey interferometric configuration, which we re- 
view in Appendix A. We also discuss spin-squeezed states 
and develop pictorial representation of those states which 
we compare to squeezed states of light and in appendix B 
we introduce the Wigner representation for a particular 
class of spin squezed states. 

In Section III, we analyze a model for spin squeezing 
based on the analogy with the optical parametric oscil- 
lator. We also seek to understand the influence of loss 
processes on the coherent spin-squeezing interaction and 
the way in which it limits the correlations achievable for a 
given interaction rate. The model consists of two bosonic 
modes (a "spin up" state and a "spin down" state) with 
loss rates and a coherent interaction that transfers pairs 
of àtoms from one mode to the other. For our simple 
model, analytical results can be obtained in the pertur- 
bative regime of small number of excitations (most àtoms 
in the lower state) and low loss rate. We estimate the con- 
ditions for which Hcisenberg-limited spin squeezed states 
can be produced. 

In Section IV, we present a scheme for inducing effec- 
tive coherent interactions between àtoms in an atòmic 
ensemble. These coherent interactions lead to massive 
entanglement of the ensemble and to characterize the 
degree of entanglement thus obtained, we calculate the 
squeezing or reduction in fluctuations of one particular 
observable. The coherent interaction is based on Ra- 
man scattering into a cavity mode for which the atòmic 
medium is made transparent by Electromagnctically In- 
duced Transparency (EIT). The slowly propagating mode 
is then best described by a polariton: a collcctive excita- 
tion that is partly photonic and partly "spin" excitation 
of the atòmic ensemble (the up and down states of the 
spin being two metastable states). The overall process 
leads to the creation of pairs of excitations, one being a 
"spin flip" created by Raman scattering, the other being 
a polariton which can be "steered" into a photon or spin 
flip excitation "on demand". We find that substantial 
spin squeezing can be obtained for atòmic ensembles in 
low finesse cavities, without the strong coupling requirc- 
ment of cavity QED. In the limit of unity finesse this 
corresponds to free-space configuration and substantial 
correlations can still be produced in this case. In the op- 
posite limit of high finesse, very strong correlations are 
obtained and in particular we estimate the regime for 
which Heisenberg limited spin-squeezed states are pro- 
duced. 
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II. RAMSEY SPECTROSCOPY WITH 
CORRELATED ÀTOMS 

In appendix A, Ramsey spectroscopy is reviewed and 
in particular we show how the phase accuracy in phase 
estimation based on the Ramsey fringe signal is, at the 
maximum sensitivity point, given by 



AJ X 



(D 



where A J x is the variance in the x-component of the 
pseudo angular momentum (of length J = N/2) repre- 
senting the state of N two-level àtoms and ( J z ) is the 
expectation value of the z-component of the pseudo an- 
gular momentum (both the variance and the expectation 
value are calculated in the initial state). 

For an uncorrelated state of àtoms e.g. with all àtoms 
in their lower state so that the state of the ensemble 
is described by \J Z = —N/2), it is found that AJ X = 

= v/l/ÏV. In order 



yJ/2 and (J z ) = —J so that 5(, 
to improve the phase accuracy, one must use a state for 
which the variance in J x is reduced while (J z ) is little 
changed. Consider therefore a state such as an eigenstate 
of J x , for example \J X — 0). Calculating the expectation 
vàlues and variances we find ( J x ) = (J y ) = (J z ) = 0, 
AJ X = and AJ y = AJ Z = ^JJ{J + l)/2. However the 
Ramsey signal has amplitude proportional to (J z ) and 
therefore vanishes for all phase angles <fi, which means 
that even though the noise or fluctuation properties of 
the signal may be improved, its average is zero. Note that 
this is because we have chosen J z ((f>) as our observable, 
other observables such as J 2 (</>) for example may lead to 
non-zero average signal together with reduced variance 
|^9|,^0|. However, it turns out their signal to noise ratio 
is very much reduced compared to that of the Ramsey 
scheme f2Cj| . It is thus necessary to consider states that 
lead to a reduced variance AJ X while maintaining a large 
signal amplitude, i.e. a large (J z ). We therefore consider 
states such as 
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where a is a real number parametrizing the state 
\tfj(a)). It is straightforward to calculate the expectation 
vàlues 

(4) = o 

(Jy) = 



(J Z ) = 
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FIG. 1. Number of àtoms detected in the upper state 
(N e ) relative to total number of àtoms (the total number is 
N = N e +N 3 and thus N e /{N e +N g ) = (N/2+{J g (<j>)))/N) vs. 
accumulated phase (f> — (u — u>o)T for a) uncorrelated àtoms 
and b) correlated àtoms in a spin-squeezed state {^(aj), for 
N = 100 and a = —1.1 (note that error bars have been mag- 
nified by a factor of 10 for clarity) . Note how squeezing of the 
variance improves the phase accuracy. 

The signal amplitude which depends on (J z ) can thus 
be rather large (0[N]) while the noise amplitude charac- 
terized by AJ X is minimized (0[1]). The Ramsey signal 
and phase accuracy for such a state is shown in Fig. lb, 
compared to the case of uncorrelated àtoms (Fig. la). 

These states are minimum uncertainty states i.e. 
(AJ x )(AJy) = h\(Jz)\ f° r all vàlues of the parameter 
a. Also, their phase accuracy is given by 



5</>(±tt/2) = 
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(5) 



which is of order Í/N. The best phase accuracy is 
obtained for a — * in which case the optimal phase ac- 
curacy is 5<j)(±7r/2) = \/2/N. Note that in this case 
the signal amplitude (oc (J z )) becomes vanishingly small 
((J z (4>)) —> f° r au 0) and also the range of vàlues of cj> 
for which improved phase accuracy is achieved becomes 
vanishingly small around <p — ±7r/2. For these reasons, 
the optimally spin squeezed state \ip(a = 0)) may prové 
impractical. Note however that for finite a i.e. for \a\ = 1, 
the signal amplitude is large (~ N/^/8) and the phase ac- 
curacy is independent of <f> 
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i.e. twice the Heisenberg limit. 
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FIG. 2. Number of àtoms detected in the upper state vs. 
accumulated phase <f> = (lú — ujo)T for correlated àtoms in var- 
ious spin-squeezed states \ip(a)), for N = 100 and a = —0.9, 
a = — 1 and a — —1.1. Also shown is the phase accuracy 
5(j>{4>) vs. accumulated phase (the dashed line represents the 
Standard quantum limit 5<j> = 1/\ÍN). Note how 5<f>(4>) gets 
to a minimum value of order 2/N = 0.02. 

In Fig. 2 we show the signal and variance for various 
spin squeezed states along with the phase accuracy 5(j){4>). 

We can gain a better understanding of the squeezing 
in the states (||) \^{a)) by looking at various representa- 
tions of them. The simplest representation is to project 
the state onto eigenstates of the three components of the 
angular momentum 



Pi (m) = \{J t =m|V(o)) | s 



(7) 



where | Ji — m) is the eigenstate of the i-component of 
angular momentum with eigenvalue m. 

From Fig. 3 it is clear that the expectation value of 
J x and J y are zero in such a state, whereas (for a = — 1) 
the expectation value of J z is large and negative, the 
variances are clearly given by (Q). It is interesting to 
note the similarity of these angular momentum squeezed 
states and those of a harmònic oscillator (i.e. squeezed 
states of light). In both cases, the probability distribu- 
tions vanish for odd number of quanta (note that for 
simplicity we consider only N even here). 
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FIG. 3. Projection of the state \ip(a)) onto eigenstates 
of the angular momentum operators J x , J y and J z for 
J = N/2 = 20 and a = -1. The P x (m) distribution is 
sharply peaked since the state \ip(a)) is a superposition of the 
m x — — 1,0, +1 components only; the P y (m) distribution is 
broad and symmetric; the P z (m) distribution vanishes for m 
odd and the even components decrease roughly exponentially 
with m. 

For even number of quanta, the behaviour is nearly 
exponential P z (m) cx e _c (" l+ ^ for some constant c. A 
mechanism for generating such states starting from the 
uncorrelated state \J Z — — J) must therefore be one in 
which àtoms are excited in pairs, i.e. two àtoms in 
the ground state are transferred to the excited state 
\9)\d) ~ * l e )l e )- Consider the similarities with squeezed 
states of light: in particular the photon number distri- 
bution vanishes for odd photon number in the case of 
squeezed vacuum due to the form of the squeezing Hamil- 
tonian H — — ix[à^ 2 — à 2 ], which creates and destroys 
photons in pairs. Since we find a similar cancellation of 
the probability of there being odd number of excitations 
for the states |^(a)), the interaction giving rise to such 
states starting from all àtoms in their lower states must 
likewise create and destroy excitations in pairs and thus 
be of the form H = —íhx[L\ — L 2 J\- This process can 
also be viewed as a coherent collision mechanism. Morc- 
over, for the whole atòmic ensemble to become entangled 
(not just particular atom pairs), this process must oceur 
completely symmetrically for all àtoms. It should not be 
two particular àtoms that get transferred to the excited 
state, rather it should be two collective excitations that 
get created 
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which is the state obtained by letting j? operate on 
\J Z = — J). The simplest Hamiltonian giving rise to this 
type of interaction is analyzed in section III to quantify 
the squeezing generated by this mechanism. This form 
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of interaction was considered by Kitegawa and Ueda || 
in their clàssic study of spin squeezing and was dubbed 
the "two axis countertwisting" interaction. 

III. TWO AXIS COUNTER TWISTING MODEL 

We now turn to thc analysis of the two axis counter- 
twisting Hamiltonian Q 

= hX^L X Ly + LyLa^j . (9) 

As argued in last section it is this type of Hamiltonian 
that most closely parallels squeezed state generation for 
light. 

We now present a general theory that allows to quan- 
tify atom correlations and takes into account decoher- 
ence and finite system size. Specifically, we consider 
two bosonic modes (such as for a two component Bosc- 
Einstein condensate or for an atòmic ensemble with two 
relevant atòmic levels) with annihilation operators ói and 
à 2 - The system is also subject to damping i.e. loss of 
àtoms at rates which may depend on the internal state. 
The equations of motion for the two modes are then (with 
L + = à\à\ and L_ = à\à 2 ) 

ài = —71 fií + xà\àl + Fx(t) 

à 2 = -72Ó2 - Xà\à\ + F 2 (t) (10) 

where Fj (t) are delta-correlated Langevin noise forces 
with appropriate diffusion coefhcients = (Fi (t)Fj (t)) . 

In order to discuss spin squeezing, it is easier to rewrite 
the equations of motion in terms of the Stokes parameters 

Lq = N = à\à\ + à\à 2 
L x = (à\ài + à\à 2 )/2 
L y = (à\à\ — à\à2)/2i 

L z = (àlà 2 -à\à 1 )/2 (11) 

for which the equations are 

Í = -2TL + A 1 L Z + F (t) 

L x = -2TL X + x(L x L z + L Z L X ) + F x (t) 

L y = -2TL y - x(L y L z + L z L y ) + Fy(t) 

L z = -2TL Z + 7L0 - 2 X (L 2 X ~L 2 y )+F z (í) (12) 

where T = (71 + 7a)/2, 7 = (71 - 72)/2 and Fj(t) 
are new delta-correlated noise forces associated with the 
damping. 

Since ( p^ ) are non-linear operator equations, in the 
equations of motion for the first order moments (Li) 



there are terms that depend on those first order moments 
but also terms depending on the second order moments 
(L{Lj). Similarly the equations of motion for the second 
order moments depend on themselves and also on the 
third order moments, and so on, leading to thc BBGKY 
hierarchy of equations of motion for the expectation vàl- 
ues of operator products. In order to solve this set of 
equations, the hierarchy must be truncated at some or- 
der pi] ]. Keeping the first and second order moments, we 
truncate thc BBGKY hierarchy by the approximation 

(LiLjLk) » (LiLj)(L k ) + (LjL k )(Li) + {t k U){tj) 

-2(L i )(L j )(L k ). (13) 

The equations of motion for the expectation vàlues 
li = (Li) and the second order moments = (LiLj + 
LjLi) — 2(Li)(Lj) are then obtained from (|Ï2]). We are 
interested in the case when all àtoms start in mode 1, 
i.e. lo(0) = NJ x (0) = l y (0) = 0,1,(0) = -N/2 and 
&xx(0) = A TO (0) = N/2 (all other second moments van- 
ish) and for simplicity we take 71 = 72 = T, 7 = 0. Writ- 
ing only the relevant equations and omitting vanishing 
terms (such as those proportional to A xz and A yz which 
are zero for all times), we have (after some àlgebra) 

lo = ~2TIq 

l'Z 2ï^l Z X(A XX Ayy) 

À xx = -4TA XX + Tlo + 4 X 1 Z A XX 

Ayy = —4TAyy + Tl - ^XlzAyy (li) 

and l x (t) = l y (t) = 0. These equations are non-linear 
and cannot be solved analytically nor perturbatively in 
T/x. For short enough times, the number of excitations 
into mode 2 is small and l z ~ —N/2, so that plugging 
this in U) we havc 

A xx (t)^^e- 2N * + 0[T/x] 

A yy (t)-^e 2N ^ + 0[r/ X ] (15) 

i.e. the variance of the x-component of the pseudo- 
angular momentum is squeezed while that of the y- 
component is anti-squeezed. Plugging these back in the 
equation of motion of l z , we obtain l z (t) ~ —N/2 + 
(cosh2Axí - l)/2 + 0[T/x\- This equation predicts 
growth of l z without bound, however we know that be- 
cause l z is the z-component of an angular momentum 
vector, we must have \l z \ < N/2. The phase space of 
this angular momentum vector is the Bloch sphere and 
in essence we have neglected the small curvature of the 
Bloch sphere (of radius R = ^V/2) and have approxi- 
mated the phase space by the flat planar phase space of 
a harmònic oscillator. We call this approximation the 
bosonic approximation, since it predicts infinite squeez- 
ing in the long time limit and in the absence of dissipa- 
tion, similar to the case of squeezed light. Formally this 
is equivalent to assuming 
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0^02 — à\à\ 
-N 



(16) 



i.e. the operator S + = à\à\/y/N obeys bosonic com- 
mutation relations. Under this approximation the Hamil- 
tonian (g) becomes H = -i{%xN/2){S\ - S 2 _) which is 
idcntical to the Hamiltonian describing squeezing of light 

In order to take into account the curvature of phase 
space and the non-bosonic nature of the angular momen- 
tum operators, we use the following transformation 



N = Nh 
L x = VNh x 
L y = VNhy 

N - 

L z = h z — —ho 



(17) 



in terms of which the commutation relations become 



x,y,z-> ^0 


= 


h z , h± 


= ±h± 


h + , h— 


hz 
~ N 



(18) 



where h± = h x ± ih y 



In the limit N — > oo these 
commutation relations become those of bosonic operators 
i.e. limAr_ +oo [/i , h z , h + , h-\ — [í,wà,à\à], a process 
formally known as a group contraction The linear 
transformation of operators ( |Ï7| ) does not introduce any 
extra approximation. 

The Hamiltonian (m can be re-expressed as 



H 



h X N (h x 
M 



hyh x 



(19) 



where we have defined £ = %N- We can now obtain 
equations of motion for the expectation vàlues hj = (hj) 
and the second order moments <5y = (hihj) — 2(hi)(hj) 



from (|14[). Letting r = N\t = & be a rescaled time, 
k = T/(Nx) = r/£ be the rescaled dissipation rate and 
writing e = í/N and x = dx/dr, we have 

ho = — 2k/io 

h z = 1k,H z (S X x àyy) 
S X x = -4k5 xx + nh - 2h S xx + Aeh z 5 xx 
5 yy = —Anòyy + nho + 2hoS yy — Aeh z 5 yy . (20) 

Note that these equations are formally equivalent to 
(|Ï4|), no approximation has been made from ( |Ï4| ) to (20). 
Letting e — > reproduces the results of the bosonic ap- 
proximation obtained above in the limit of l z — —N/2. 
Terms of order e and higher represent corrections to the 



bosonic approximation and, as shown below, they give 
rise to a limit to the amount of squeezing achievable. 

Solving ( |2C| ) to first order in e and k we obtain, writing 
only the relevant terms, 
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Mr) = - [e~ 2T + ( K + e/2) + 



(21) 



which shows that the variance AJ X = d (N/2)5 XX is 
squeezed. Second order terms in k and e come multi- 
plied by an exponentially growing term e 2r so that as a 
function of time, the variance reaches a minimum valuc 
S xx ~ max[/í, e] at a time e _T * ~ max[ft, e], after which it 
grows exponentially and the squeezing is lost. Note that 
this behaviour (5 XX (t) reaches a minimum value and then 
increases again) also occurs when k — > 0, indicating that 
it is a genèric feature of the finite system size. This model 
predicts that a variance S xx ~ e = Í/N — * AJ X ~ 1 is 
achievable as long as losses are small enough, i.e. k ~ e, 
which in terms of \ and T means 



X ~ T or C ~ NT 



(22) 



where x corresponds to the single-atom nonlinear in- 
teraction rate and T represents the single-atom loss rate. 

In order to achieve any squeezing (S xx < í/2) it is 
necessary to have n < 1 i.e. Nx > T. In the regime 
Nx 3> r very strong correlations can be obtained. Note 
that the single-atom nonlincarity can still be relatively 
weak compared to the single particle loss rate (x <S T). 
For example when the dissipation rate is such that k ~ 
yfk i.e. VNx ~ r, the amount of squeezing obtained ( pï| ) 
is S xx ~ í/\/N. It takes a time e~ 2r ~ \fi to reach this 
state and the number of partides lost during that time is 
AiV ~ N x 2kt ~ v^VlogA^. This number can therefore 
also be thought of as the maximum number of partides 
that can be lost from the ensemble without destroying 
squeezing beyond S xx ~ 1/ y/N. 

In order to reach the Heisenberg limit it is required that 
the single atom nonlinearity x be larger than the decay 
rate T. Note that in this case, the number of àtoms lost 
by the time optimal squeezing is achieved is AN ~ log N 
which indicates that a very small number of àtoms is lost. 
This number also corresponds to the maximum number 
of partides that can be lost from the ensemble and not 
destroy squeezing at the Heisenberg limit level. Clearly 
the more squeezed the state of the àtoms is, the more 
sensitive it becomes to atom loss and in general to any 
form of dissipation. 



IV. COHERENT ATOM INTERACTIONS VIA 
SLOW LIGHT 

We now describe a technique to induce effective coher- 
ent interactions between àtoms in metastable states Hi. 
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The technique is based on a resonantly enhanced nonlin- 
ear process involving Raman scattering into a "slow" op- 
tical mode ]Í5| , which creates a pair of spin-flipped atom 
and slowly propagating coupled excitation of light and 
matter (dark-state polariton). When the group velocity 
of the polariton is reduced to zero Jï(],[ïj]] , this results in 
pairs of spin fiipped àtoms. The fact that pairs of atòmic 
excitations are created in this process can also be viewed 
as a coherent interaction between àtoms, i.e. a controllcd 
"collision" lcading to entanglement of the state of each 
atom with that of every other atom in the ensemble. 

A number of proposals have been made for gener- 
ating entangled states of atòmic ensembles and result- 
ing in so-called spin squeezed states. Some are based 
on interatomic interactions at ultra-cold temperatures 
p3|, whereas others involve mapping the states of non- 
classical light fields into àtoms 0] , QND measurements 
of spins |25|| with light or dipole blockade for Rydberg 
àtoms [|26f. In contrast to some of these mechanisms 
the present approach does not require coherence of the 
atòmic motion or sources of non-classical light and is 
completely deterministic thereby significantly simplifying 
possible experimental realizations. We further show that 
the present technique can be made robust with respect 
to realistic decoherence processes such as spontaneous 
emission and leakage of slow photons from the medium. 

We considcr a system of N àtoms (Fig. 4a) interacting 
with two classical driving fields fi^a and one quantizcd 
mode à of a running wave cavity that is initially in a vac- 
uum state. Note that we consider a cavity configuration 
for ease of theorerical treatment, the results of this anal- 
ysis however remain vàlid in the limit of unity finesse, i.e. 
in free space configuration. Relevant atòmic sublevels in- 
clude two manifolds of metastable states (e.g hyperfine 
sublevels of electrònic ground state) and excited states 
that may be accessed by optical transitions. The àtoms 
are initially prepared in their ground states |<7). One 
of the classical fields, of Rabi frequency Oi, is detuned 
from the atòmic resonance by an amount roughly equal to 
the frequency splitting between ground state manifolds. 
The other field of Rabi frequency ÇI2 is resonant with an 
atòmic transition | £>2 ) — > {0-2) ■ The quantized field can 
be involved in two Raman transitions corresponding to 
Stokes and anti-Stokes processes. Whereas the former 
corresponds to the usual Stokes scattering in the forward 
direction, the latter establishes an Electromagnetically 
Induced Transparency (EIT) and its group velocity is 
therefore substantially reduced. 




FIG. 4. Energy levels scheme for the efïective coherent 
interaction leading to creation of pairs of àtoms a) in differ- 
ent final states ( "non-degenerate" scheme) and b) in identical 
final states( "degenerate" version) . 

The pair excitation can be viewed as resulting from 
quantized photon exchange between àtoms (Fig. 5) in 
a two-step process. The first fiipped spin is created due 
to Stokes Raman scattering, which also results in photon 
emission in a corresponding Stokes mode. In the presence 
of EIT, this photon is directly converted into a dark-state 
polariton which becomes purely atòmic when the group 
velocity is reduced to zero. This implies that atòmic spins 
are always fiipped in pairs. In Fig. 4a the two final states 
involved in Raman transitions are different and atòmic 
pairs in different states are created. In Fig. 4b the final 
states of the two Raman processes are identical, in which 
case atòmic pairs in the same state result. The analysis 
of this "degenerate" version of the scheme is similar to 
the "non-degenerate" case and we will consider only the 
latter case here. 




FIG. 5. Diagram illustrating coherent atom-atom interac- 
tion mediated by dark-state polariton, leading to the creation 
of a pair of spin-flipped àtoms. 

For conceptual simplicity we assume that the quan- 
tized field corresponds to a single mode of a running- 
wave cavity with a creation operator and atom-field 
coupling constants 51 and g 2 - The interaction Hamilto- 
nian for the system of N àtoms and light can be split into 
two parts H = Hji am + H res corresponding to the Stokes 
Raman process and the anti-Stokes process respectively: 

HRam = —Tl&^alal ~ ^iT,b l0l 

+ [nn^g + ngiaz aibl + h.c], (23) 

H res = hS 2 T,b 2 b 2 + T%52^a 2 a 2 

+ [% 2 as a29 + nn 2 j: a2b2 + h.c], (24) 

where = |/lí)íí(i/| are collective atòmic opera- 
tors corresponding to transitions between atòmic states 
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\i/), A is the detuning of the classical field íii from 
the single-photon transition \g) — > \a%), 8\ and 8 2 are the 
two-photon detunings from the \g) — > |6i) and |<?) — > |ò 2 ) 
transitions respectively as shown in Fig. 4. 

In the limit of large detuning A and ignoring two- 
photon detunings for the moment, the Hamiltonian 
HR am describes off-resonant Raman scattering. We take 
into account realistic decoherence mechanisms such as 
spontaneous emission from the excited states in all di- 
rections and decay of the cavity mode with a rate k. 
The evolution of atòmic operators is then described by 
Heisenberg-Langevin equations: 



H, 



eff 



-iH{P D S\ - SíPd), 



(29) 



(25) 



where 7^ is a decay rate of coherence [i — > v and 
are associated noise forces. The latter have zero average 
and are <5-correlated with associated diffusion coefficients 
that can be found using the Einstein relations. 

After a canonical transformation corresponding to adi- 
abatic elimination of the excited state (see Appendix C 
for details), Hn am becomes equivalent to the effective 
Hamiltonian 



HRam = %xà^Sl + h.c. 



(26) 



where Si = Y, gb J^N and X = gWNÜl/A. This 
effective hamiltonian thus describes the process in which 
a Stokes photon is emitted necessarily accompanied by a 
spin flip. The quantum state of the Stokes mode is thus 
perfectly correlated with the state of the atòmic spin flip 
mode. 

The resonant part of the Hamiltonian H res is best an- 
alyzed in terms of dark and bright-state polaritons [E7| 



Pd 



Pb = 



ü 2 a - g2vNS 2 
g 2 VNa -i 



-01 



•«I 



(27) 



which are superpositions of photonic and atòmic excita- 
tions à and S2 = ^ g b 2 / V~N.. In particular, H res has an 
important family of dark-states: 



\D r< 



(4n<?)|vac) 



(28) 



with zero eigenenergies. This means that once in the dark 
state, the system stays in the dark state. Note that all 
other eigenstates of H res have, in general, non-vanishing 
interaction encrgy. Under conditions of Raman resonancc 
and sumciently slow excitation ( "adiabatic condition' , 
see Appendix D for details) the Stokes photons emitted 
by Raman scattering, Eq. ( |26| ) , will therefore couple solely 
to the dark-states (28). In this case the coherent part 
of the evolution of the entire system is described by an 
effective Hamiltonian: 



with £ = QiQ 2 /A x giVN '/y/g$N + flj (without loss 
of generality, £ was chosen imaginary here for simplified 
calculations) . The Hamiltonian ( p9j ) describes the coher- 
ent process of generation of pairs of excitations involving 
polaritons and spin-flipped àtoms. Note that for small 
numbcr of excitations the spin waves and polaritons obey 
bosonic commutation relations and this Hamiltonian is 
formally equivalent to that describing optical parametric 
amplification (OPA) of two modes In the non-bosonic 
limit, this Hamiltonian is also analogous to the "counter- 
twisting" model of (||). In appendix D we show that the 
coupled equations for the polariton Pd and the spin flip 
S'í are given by 



Sí = (g^7L - 1L - i8x)S\ + ÍP D + Fl (t) (30) 
P D = -{k/t, + ió 2 )P d + £SÍ + F D {t) (31) 

where the polariton decay rate includes an atòmic part 
7l and a photonic part k/ï) due to leakage of photons 
out of the medium (at a rate reduced by the factor 
V = | <?2 1 2 ^V/ 1 02 1 2 equal to the ratio of vacuum light ve- 
locity to the group velocity of slowly propagating Stokes 
photons). The spin flip operator equation ( |30| ) is seen to 
contain both a decay term and a gain term due to sponta- 
neous emission into the bright polariton mode. Note that 
this apparent decrease in dissipation is however accompa- 
nied by increased fluctuations denoted by the new noise 
force operator Fg^í). The effective detuning between 
the polariton and spin flip mode is seen to correspond to 
the difference in two-photon detunings 61 —62. 

We now consider the scenario in which the system is 
evolving for a time r under the Hamiltonian H e f f , after 
which both fields are turned off. If the procedure is adia- 
batic upon turn-off of the coupling fields fli t2 the polari- 
tons are converted into pure spin excitations Pd — > S 2 . 
Hence the entire procedure will correspond to the follow- 
ing state of the system: 

I*) = ^>7E( tanh ^) n ^(^) n ( 5 í) ,l |5)|vac) 

n 

^ ^77Z)( tanh ^ r ) n | n 6i- n fe)|vac). (32) 

n 

Here \n bl ,n b2 ) = l/n\{S^) n {S^) n \g) are Dicke-like 
symmetric states of atòmic ensemble and we assumed 
n biM ^ N. For non-zero £t this state describes an en- 
tangled state, for which relative fluctuations between the 
two modes decreases exponentially to vàlues well below 
the Standard quantum limit (SQL) corresponding to un- 
correlated àtoms. 

The present technique can also be viewed as a new 



ït 



mechanism for coherent "col·lisions" |13| between àtoms 
mediated by light. In particular, the case when atòmic 
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pairs are excited into two different levels (as e.g. in Fig. 
4a) closely resembles coherent spin-changing interactions 
that occur in degenerate atòmic samples S , whereas the 
case when atòmic pairs are stimulated into the identical 
state (Fig. 4b) is reminiscent of dissociation of a molec- 
ular condensate |29| . To put this analogy in perspective 
we note that the rate of the present optically induced 
process can exceed that of weak interatomic interactions 
by orders of magnitude. Therefore the present work may 
open up interesting new possibilities for studying many- 
body phenomena of strongly interacting àtoms. 
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FIG. 6. (a) Quadrature variance AY+ 



single-photon 



detuning A and interaction time £í, (b) same for A = A opt 
and Si = 82 showing maximum squeezing AY+ ~ 0.02 (for 
yj g 2 N /jk = 100), (c) Number of excitations pumped in the 
system vs. time (same conditions as in b) and (d) AY+(i*) 2 
vs. two-photon detuning 8 = (81 — 82) /2 for A = A op t and 
where í* gives maximum squeezing. 

To quantify the resulting correlations established be- 
tween the polariton mode Pr> and the pure spin flip mode 
S'í , we introduce the quadratures of both modes (which 
are bosonic for small number of excitations) in direct 
analogy to the optical parametric case. We define the 
quadratures X 1 = {S x + Sf)/V2, Y\ = i{S\ - Sf)/V2 
and similarly for the polarition mode; these can be mea- 
sured e.g. by converting spin excitations to light. Cor- 
relations between the modes appear due to dynamical 
evolution and squeezing is found in certain quadratures 
of the sum and difference modes A"_ = (Xi — Xd)/V% 
and Y + = (Y\ + Yd)V2- In the language of harmònic 
oscillators, the positions in mode 1 and 2 = D are corre- 
lated (X\ ~ Xd), while the momenta are anti-correlated 
(Yi ~ — Yd). For small number of excitations the sum 
and difference modes obey Standard commutation rela- 
tions [X Q ,Yg] — ~iS a ,f3 where a,(3 = +, — or 1,13. A 
quadrature Y± is squeezed when AY±(í) 2 < 1/2 and the 
Heisenberg limit corresponds to AY±(t) 2 ~ 1/N. 

We find that squeezing is optimal under conditions 
of four-photon resonance (<5i = 62) and in the limit of 
77 ^ 1 (Fig. 6). Evolution leads to squeezing of Y+ and 
AT_ , anti-squeezing of Y_ and X + . The squeezing in Y+ 



reaches a minimum value at t = t* after which the grow- 
ing fluctuations in X + give rise to increased noise in Y + . 
Note that the total number of excitations (both modes) 



in the system, equal to (Xj_ + X_ 



Y_), grows ex- 



poncntially with time (Fig. 6c). Specifically, in the case 
.9i = 92 = 9 and thus 71 = 72 = 7, for Çt > 1, we have: 



(AY + (Í)) 2 = 1/2 {e- 2 ^ - 
' 1L + k/v 



37l + k/v 



(33) 



2Í 



where we have neglected terms of higher order in 7^ /£ 
and k/(t)Ç). The maximum amount of squeezing is ob- 
tained after an interaction time t* such that e~ 2 ^* = 
{lL + k/t])/2^ and is given by 



(AY+) 2 = 



5j L + 3k/t] 
4£ 



(34) 



i.e. of the order of the damping rate divided by the 
coherent interaction rate. 

Since both the interaction parameter Ç and the relax- 
ation rate of the polariton jd = 7l + t/77 depend on 
the single photon detuning A (Fig. 6a), we find that 
squeezing is optimized for 



A opt = 71 



/ 5|»!| 2 |g| 2 iV 
3|íí 2 | 2 7« 



(35) 



and with this optimal valuc of the detuning, the squeez- 
ing reaches a minimum value of 



y/1574 
V\9\ 2 N/ik 



(36) 



Note that the denominator is equal to the atòmic density- 
length product multiplied by the empty cavity finesse 
and can easily exceed 10 4 even for modest vàlues of the 
density-length product and cavity finesse. We further 
emphasize that typical generation rate resulting in such 
optimal squeezing Q1Q2/ A opt can easily be on the order 
of fraction of MHz. In such a case other decoherence 
mechanisms are negligible. Doppler shifts can also be 
disregarded as long as all fields are co-propagating. 

For the "degenerate" version of the interaction (i.e. 
with identical final states for the spin fiips, see Fig. 4b), 
the effective hamiltonian can be written as 



H eff =ihÇ(SV-S 2 ) 



(37) 



where the limit 77 ^ 1 has been used to write Pjj ~ —5, 
with S = 1/VÏVS 9 6 the spin flip operator. In this case 
the correlations lead to squeezing of X = (S + S*)/y/2 
and anti-squeezing of Y = i{S-S^)/y/2. The analysis for 
this configuration is very similar to the non-degenerate 
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version, in particular the maximum amount of squeezing 
achievable is also given by an expression of the form (36). 

We can now obtain a condition for achieving 
Heisenberg-limited spin squeezed states, i.e. (AY +opt ) 2 ~ 
1/N. We see from (pj) that this requires 



£ - nt 



(38) 



where T = (57l + 4k/?7)/4 is the effective damping rate 
of the system. This is in complete agreement with the 
estimate based on our simple bosonic model of section 
III (|2^). In terms of the single photon Rabi frequency g, 
the cavity decay rate k, the spontaneous emission rate 7 
and the number of àtoms N, the condition for achieving 
some squeezing i.e. (AY~ + ) 2 < 1/2 is 



\g\ 2 N>K 1 



(39) 



which can be easily achieved in the laboratory since 
it simply corresponds to the condition that the density 
length product multiplied by the cavity finesse be larger 
than one. In the cavity QED regime of strong cou- 
pling \g\ 2 ~ K7, very strong quantum correlations i.e. 
(AY" + ) 2 ~ 1/yN between àtoms can be produced. In or- 
der to obtain Heisenberg limited spin squeezed states i.e. 
(Ay + ) 2 ~ 1/iV, one requires a more stringent condition 



1.91 



Nn-y 



(40) 



which can be fulfilled only in the strong coupling 
regime of cavity QED for a limited number of àtoms. 
Note that this regime has been achieved experimentally 
by several groups [ pÜ| and would allow for Heisenberg 
limited spin squeezing for as many as ~ 10 3 àtoms. We 
have analyzed in this paper the situation of a running- 
wave cavity, so that all àtoms couple equally apart from 
a possible phase to the cavity mode irrespective of their 
position. In order to fulfill the cavity QED regime, small 
cavity volume is needed i.e. standing wave cavities. For 
àtoms in such a cavity the coupling to the cavity mode 
is position dependent and it becomes necessary to lo- 
calize àtoms accurately at the antinodes of a trapping 
mode. Note that significant experimental progress has 
been made towards this direction by several groups [^ï] . 
Once the àtoms are well localized in the cavity, the in- 
teraction can proceed via a neighbouring mode b (e.g. 
different from the trapping mode a) so that for àtoms lo- 
calized within a small region in the cavity the two modes 
have essentially the same wavelength and àtoms would 
therefore couple equally to the b mode as well, irrespec- 
tive of their position. 



V. DISCUSSION AND CONCLUSION 

We have reviewed Ramsey spectroscopy and the use of 
spin squeezed states in precision measurements of this 



type. With the experimental motivation of minimiz- 
ing the phase accuracy in phase estimation with Ram- 
sey fringes, we introduced a particular class of squeezed 
states. These states lead to Heisenberg limited phase ac- 
curacy and we developed various pictorial representations 
for them. The strong similarities of these representations 
of spin squeezed states to those of squeezed states of light 
suggests an analogy extending to the type of interaction 
that gives rise to squeezing. We are thus lead to consider 
the so-called "counter twisting" Hamiltonian, which has 
been shown to lead to maximal spin squeezing. We have 
studied this model for spin squeezing in the presence of a 
dissipation mechanism and analyzed the effect of damp- 
ing and finite system size on the amount of squeezing 
achievable with such an interaction. The analysis was 
based on a decorrelation approximation to the BBGKY 
hierarchy of equations of motion, followed by the use of a 
linear transformation which in the limit of large number 
of àtoms 1/N — > "contracts" the angular momentum 
operators onto bosonic operators. This allows for the 
systematic inclusion of finite system size effects. It ap- 
pears that Heisenberg limited spin squeezed states may 
be produced when the single atom nonlinearity exceeds 
the single atom loss rate. In this case the maximum 
number of àtoms that can be lost before quantum correla- 
tions are destroyed to the point of compromising the spin 
squeezing is of the order AN ~ log N. For spin squeez- 
ing at a more modest level than the Heisenberg limit, 
larger number of àtoms may be lost without compromis- 
ing the squeezing, indicating the stronger sensitivity of 
spin squeezed states to dissipation for larger amounts of 
squeezing. 

We have also presented in detail a scheme based on the 
interaction of coherent classical light with an optically 
dense ensemble of àtoms that leads to an effective coher- 
ent spin-changing interaction involving pairs of àtoms. 
Àtoms may be transferred to the same final state lead- 
ing to spin squeezing (analogous to squeezing of light by 
degenerate OPO) or to different final states in this case 
leading to quantum correlations between different atòmic 
modes (analogous to quantum correlations between elec- 
tromagnètic modes by non-degenerate OPO). We have 
shown that this process is robust with respect to realistic 
decoherence mechanisms and can result in ràpid genera- 
tion of correlated (spin squeezed) atòmic ensembles. The 
amount of correlations created by this effective interac- 
tion can be simply expressed in terms of the single photon 
Rabi frequency <?, the atòmic spontaneous emission rate 
7 and the cavity decay rate k. We find that the gener- 
ation of spin squeezed states requires g 2 N ~ «7, which 
can easily be achieved in low finesse cavities with e.g. 
room temperature atòmic vapours. Very strongly corre- 
lated states can be produced when the strong coupling 
regime g 2 ~ kj of cavity QED is achieved and the gener- 
ation of Heisenberg limited spin squeezed states requires 
g 2 ~ Nivy. The effective interaction rate £ = Íiií^/A 
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which depends on the Rabi-frequency of two applied clas- 
sical fields and a detuning from an atòmic transition 
A can be fast and is controllable. Furthermore, the re- 
sulting spin excitations can be easily converted into pho- 
tons on demand, which facilitates applications in quan- 
tum information processing. Possible applications involv- 
ing high-prccision measurements in atòmic clocks can be 
also foreseen. 

We thank M.Fleischhauer, J.I.Cirac, V.Vuletic, S.Yelin 
and P.Zollcr for helpful discussions. This work was sup- 
ported by the NSF through the grant to the ITAMP. 



APPENDIX A: RAMSEY SPECTROSCOPY 

In Ramsey spectroscopy [0 , a collection of N two-level 
àtoms are madc to intcract with two separated fields (in 
time or in space). The lower and upper states (refered 
to as ground and excited state) have an energy diffcr- 
ence Tiloq and àtoms will thus acquire a different dynam- 
ical phase e~ lEt / h depending on which state they are in. 
The effect of properly chosen electromagnètic fields is to 
perform a transformation that prepares the àtoms in a 
superposition of the two states \g) and |e). The differ- 
ent parts of the wavefunction of àtoms (corresponding to 
the ground and excited state) acquire a relative phase 
due to dynamical evolution and when the inverse trans- 
formation is applied, an interference effect is obtained. 
An exact parallel with the Mach-Zendcr intcrfcrometer 
can be drawn p^ |: the transformation preparing àtoms 
in a superposition of ground and excited states is equiv- 
alent to the transformation that lets a photon incident 
on a beam splitter explore the two arms of an interfer- 
ometer. The relative phase acquired in the two atòmic 
states during free evolution of duration T is the equiva- 
lent of the relative phase acquired by photons travelling 
in the arms of the interferometer. Finally, the second 
pulsc that performs the inverse transformation on àtoms 
is the equivalent of the recombination of signals from the 
two interferometer arms on a beam splitter. At the end 
of this sequence, the number of àtoms in either states, 
equivalent to the number of photons from either output 
of the final beam splitter, is measured. In this way, the 
signal measured depends on the acquired relative phase 
which can thus be estimated with some accuracy. 

We will now quantify this more precisely: let the fre- 
quency of the applied electromagnètic pulses be w, and 
the time delay between the two zones of interaction be 
T. The duration and strength of the applied fields are 
chosen so as to lead to tt/2 pulses, i.e. transformation of 
the atòmic state according to 



\e) 

\g) 



|e)-*|g) 
V2 

\±±j\g) 

V2 ' 



During their free evolution between the two zones, 
àtoms in the ground and excited states acquire a relative 
phase (fi which, in a frame rotating with the frequency of 
the applied field, is <fi = (lu — lo )T. 

Before entering the first interaction zone, the àtoms 
are prepared in their lower state \g) and at the exit of 
the second zone, the number of àtoms in states |e) and 
\g) is measured. 

For simplicity, we consider the case when the first 
zone leads to a tt/2 pulse and the second one a —tt/2 
pulse. The picture of angular momentum is particularly 
well suited to discuss the Ramsey interferometric scheme 
and leads to an intuitive pictorial representation of the 
scheme. The Schwinger angular momentum operators 
are defined as 



Jx = Ç^eg + ^ge)/2 
Jy = (^eg — ^ge)/2ï 

Jz = (S Ee — Ylgg)/2 



(A2) 



where £ M „ = 2~2j=i \l· i )jj( v \ are collective operators. In 



terms of these, a single tt/2 pulse ( |Al| ) is represented by a 
rotation of the pseudo angular momentum vector around 
the x-axis by an angle tt/2. For a single atom we have 
the correspondence | f) = |e) and | J.) = \g). Under a tt/2 
rotation about the x-axis, the state |) transf orms to | J y = 
-1/2) = (| T) - i\ i))A/2~ as indicated in (O. For N 



àtoms, we can think of the N individual spin ^ partides 
combining to form a pseudo angular momentum vector of 
length J = N/2. The state of the collection of N àtoms 
can then be represented by appropriate superpositions of 
the states |J, M) where — J < M < J. Of course, only 
states within the completely symmetric subspace of the 
full 2^-dimensional Hilbert space can be represented in 
this way, which is justified since the coherent interaction 
of the electromagnètic fields with the àtoms couple only 
to this symmetric subspace (i.e. all àtoms couple equally 
to the fields). 

Free evolution in the rotating frame corrcsponds to ro- 
tation of the angular momentum around the z-axis at an 
angular velocity lo — lüq. The whole Ramsey scheme can 
then be represented by the sequence: 7r/2 rotation about 
x-axis, <fi rotation about the z-axis and —tt/2 rotation 
about the x-axis. This is the transformation perfomed 
by the unitary operator 



ü(4>) = e W2J* e -^J, e -W2j* 



(A3) 



where <fi = (u> — lüq)T as before. At the end of the 
scheme, the number of àtoms in states |e) and \g) is mea- 
sured, or equivalently their difference J z (<fi) where 



(Al) 



= J z cos (f> — J x sin i 
The Ramsey signal is thus 



(A4) 
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{J z (4>)) = (J z ) cos<?f> - (J x ) sin0 
and its variance AJ z {<p) is 
AJ Z (» = [(AJ Z ) 2 cos 2 (j) + (AJ X ) 2 sin 2 

- cos <j> sin c/>({j x J x + J Z J X ) - 2{J Z ){J X )) 



(A5) 



1/2 



(A6) 



where the variance is defined as (AA) 2 = (À 2 ) — (À) 2 . 
From the signal one wants to estimate the phase <j) and 
thus the frequency difference uj — ujq. The phase accuracy 
achievable from such a measurement is related to the 
signal variance (the "noise" ) by 



AJM 

\d(JzM)\ 



(A7) 



For states such that (J x ) = (all the states we will 
consider in this paper are of this type), the sensitivity 
\d( J z ((f>)) / d(j>\ is maximal for <fr = ±ir/2 and the phase 
accuracy can be expressed as 



<^(±tt/2) = 



AJ X 

l(4)f 



(A8) 



Since AJ X and ( J z ) depend on the initial state, we see 
that different initial states lead to different phase accura- 
cies. Of particular importance is the accuracy achievable 
when all àtoms are prepared in the same initial state. In 
this case the state of the atòmic ensemble is a pure state, 
but it is however an uncorrelated state of the atòmic en- 
semble (i.e. it can be factorizcd |W) = njLi \4 ) )j)· 

Consider the case of uncorrelated àtoms for which all 
àtoms have been prepared in the lower state \g), some- 
times called a Bloch state. The state of the atòmic en- 
semble can thus be expressed in terms of eigenstates of 
the collective angular momentum operators as 



N 



Jl\g) j = \J = N/2,J x = -N/2) 



(A9) 



where J = N/2 since there are N 2-level àtoms, equiv- 
alent to N spin i partides. For such a state, the ex- 
pectation value of the angular momentum operators and 
their variances are calculated to be (J x ) = (J y ) = 0, 
(J z ) = -J, AJ X = AJ y = y/j/2 and AJ Z = 0. The 
signal and its variance are thus 



{J z (4>)) = ~J COS0 
AJ z (0) = ^/J/2sm(j). 

The maximum sensitivity is achieved at 

1 1 

7È 



^(±tt/2) = -= = — = 



2J 



(A10) 
±tt/2 

(All) 



which is the Standard quantum limit (SQL). Perform- 
ing the experiment on N independent àtoms all prepared 



in the same initial state is thus equivalent to repeating 
the experiment on one atom N times and leads to an ex- 
pected 1 / \/N factor of improvement in accuracy over the 
one atom result AS X /(S Z ) = 1. This is the best accuracy 
achievable with àtoms all prepared in the same initial 
pure quantum state. The number of àtoms detected in 
the upper state, given by {N+(<j>)) — N/2 + (J z (</))), and 
its variance are shown in Fig.la. 

There is a lowcr bound on the phase accuracy, 
set by Heisenberg's uncertainty principle, AJjAJj > 
i|([Ji, Jj])| where i,j — x,y,z. It is straightforward to 
show that 



> 



1 

N 



(A12) 



which is known as the Heisenbcrg limit. 

We now see from ( |AS| ) that in order to surpass the 
SQL, the atòmic ensemble must be prepared in a state 
such that AJ X /\(J Z )\ < l/y/N, which is a necessary 
and sufficient condition for entanglement of an atòmic 
ensemble [ p3[ . It is thus important to have a state 
for which the variance AJ X is reduced compared to its 
value for the uncorrelated state ( |A9| ) while maintaining 
a large value for ( J z ) so that the amplitude of the signal 
(J z (</>)) = (J z ) cos0 is not compromised p0| |. Such states 
which have reduced uncertainty in one observable AJ X 
(at the expense of the conjugate observable AJ y having 
increased fiuctuations) have been called spin-squcczcd 
states H. 



APPENDIX B: SPIN SQUEEZED STATES - 
WIGNER FUNCTION REPRESENTATION 



We now consider the Wigner function representation 
of the states \tp(a)). The Wigner distribution of general 
angular-momentum states [|33f is obtained from an ex- 
pansion of the density operator in terms of the multipole 
operators 



2J +k 
k=0 q=-k 

where the multipole operators are 

t kq = v v {-ly-^vw+ï ( J k J , 

y jL^í jL^í V —777, q 777/ 

7n— — J rn' — — J 



(Bl) 



and 



x \ J, m){J, m'\ 
J k J 



(B2) 

is the usual Wigner 3j symbol. 



-m q m 

The wigner distribution is thcn given by 

2.7 +k 
k=0 q=-k 



(B3) 
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where p kq = (f kq ) = Tr [pf kq ] and Yg(d,4>) are the 
spherical harmonies. In Fig. 7, the Wigner function for 
the state \ip(—l)) clearly shows the way in which this 
state has a large negative expectation value for J z , re- 
duced variance in J x and increased variance in J y . 




FIG. 7. Wigner function representation of the state \tp(aj) 
with a = —1). Plotted is the surface r(8, <j>) — W(9, <j>), show- 
ing the large and negative value of (Jz), reduced variance AJ X 
and correspondingly increased variance AJ y . 



APPENDIX C: ADIABATIC ELIMINATION OF 
EXCITED STATE IN RAMAN SCATTERING 

From the Hamitonian ( p3| ) , we obtain the equations of 
motion for the cavity mode and the ground state coher- 
ence 

à = -na - iglY, biai - ig* 2 H ga2 + F a (t) (Cl) 
E ff bi = -(7o - iSi)^ gbl + iOiE aii)1 - ïg^a t S gai 

+ F gbl (t) (C2) 

and the optical polarizations associated with Stokes 
emission evolve according to 

Éfc iai = -[7 ^ *(A - 5i)]S6 1Q1 - ïfiiSfc ig 

- i.gia(E hlfcl - S 01Q1 ) + F biai (t) (C3) 
S Q1S = -(7 + iA)S aig - ifl*[(E aiai - E 53 ) 

+ ig* lt Jv big + F aig (t) (C4) 

where we assume that population in the excited state 
|ai) decays towards |&i) at a rate 71, towards \g) at a rate 
72 and we assume a dephasing rate 70 for ground state 
coherences (7 = (71 + 72V2 and 7 ^> 70). 

We proceed by adiabatic elimination of optical polar- 
izations associated with Stokes emission. To this end we 
assume large single-photon detuning A ^> 7 and to first 
order in à we obtain (E gg ~ N) 



v - ■ 7 \ v 

^biai — l ~^) 2 -'b 1 g 

V aig = ^N(l+il)- 



. F biai (t) 

* A 

; Fq ig (t) 

A 



(C5) 
(C6) 



which we substitute in (C2) and obtain for the ground 
state spin flip operator Si = Yi gbl /^f~N 



Si = 



[(70 + 7l) i(6x + S L )} 5 X - » glV f"V 



F Sl (t) 



(C7) 



where 7/, = 7|í7i| 2 /A 2 is an optical pumping rate, 
§l = |íii| 2 /A is the light shift and Fs 1 (t) is a modified 
noise force. Light shifts can be incorporated in a redefini- 
tion of the energies and we ignore them in the remainder 
of this paper. Since the ground state decoherence rate is 
typically very small, we also assume 71, 3> 70 and in that 
limit the new <5-correlatcd noise forces have correlations 



(F Sl (t)Fl(t'))=2 7L ^S(t-t') 
7 

(Fl(t)F Sl (t'))=2 7L ^S(t-t') 

1 



(C8) 
(C9) 



APPENDIX D: ADIABATIC ELIMINATION OF 
BRIGHT POLARITON 

After adiabatic elimination of the excited state |ai), 
the relevant equations of motion are 



à 
S 2 

^ga 2 = -(7 + i5 2 )Y, ga2 - S 2 - ig 2 Na 



-(jL+iS^St + ixa + F^t) 
-na - ix*S{ - ig 2 ^ ga2 + F a (t) 
-ilL + iS 2 )S 2 - i-^E 9í2 + F S2 (t) 



+ F ga2 (t) 



(Dl) 



where S 2 = Y, gb jVN. 

From (Dl) and ( p7| ) and in the limit of large ratio of 
speed of light in vacuum to group velocity of Stokes pho- 
tons i] = |<72| 2 -^V/|^2| 2 ^> 1, we obtain the equations of 
motion in terms of bright and dark polaritons 



sl 



'D 



-(7L + ííi)Sj + i X (^- + P B ) + Fl(t) 



X* ot K ~ r 2 



Pe 



->ga 2 



-{ K / v + T 2 )P D -i^-Sl 
F D (t) _ 

-(« + t 2 /v)Pb - i X * Sil - ig 2 VNZ s 

F B (t) 



Pb 



= -rE ga2 - ig 2 VNP B + F ga2 {t) 



(D2) 

(D3) 

(D4) 
(D5) 



13 



wherer 2 = Jl+íSï, T = 7+ÍÍ2 and E sa2 = S ga2 /\/iV 
and noise forces were modified appropriately. Note that 
in the picture of dark and bright polaritons, only the 
bright polariton is coupled to the excited state through 
the optical coherence T> ga2 . 

Under adiabatic conditions, the bright polariton 
evolves slowly (on a typical timescale T) and we can 
solve perturbatively in 1/T. The equations (D4) and 
( p5| ) are of the form x = M.x + y, where x is the vec- 
tor (P B , S g a 2 )j M is a 2 x 2 matrix and y is a source 
term 



d_ 
dt 



Pb 



k ig 2 VN \ P B 
ig-2"/N r / [ Sga 2 

iX*St-JP D +F B (t) 

Fga 2 (t) 



(D6) 



where we have used n 3> 7l and where F B (t) and 
F g a 2 (t) are appropriate noise forces. Thcsc equations 



can be solved easily to first order by y^°\t) 
higher order approximations yielding x*-™) (t) — 

We can rewrite 



|g 2 | 2 ^v 

«7 



3tt x 



—LA 2 
V 



x T 



M _1 .[y— 



(D7) 



i.e. the dcnsity lcngth product multiplied by the cav- 
ity finesse, so that with densities corresponding to room 
temperature atòmic vapours, optical wavelengths and fi- 
nesse of order 100 this quantity is already of order ~ 10 4 . 
We can thus assume that \g2\ 2 N/(nj) » 1 and solve in 
powers of nj/(\g2\ 2 N). 

We see from ( |D6| ) that xW(í) is of order 
[n'y / (\g2\ 2 N)]( n+1 ï and thus solving to lowest order we 
find 



Pb = 



1 



152 



7 -N l 



-iT X *St - ^P D 



TF B {t) 



~ Íg2Fga 2 (*)] 

so that when 77 ^> 1, 

ac^ + P B 
V 



(D8) 



77 \g 2 \ 2 N l 

Íg2F ga2 (t)}. 



-iT X *St + TF B (t) 



(D9) 



The co uple d equations of mo tion for the dark state 
polariton (D2) and the spin flip (D2) then become 



Sl = (g^7L -1L- i5i)S\ + i-±=P D + Fl (t) (D10) 
P D = -{n/r, + lL + i5 2 )P D - i^-S\ + F D (t) (Dll) 



\9i\ 



■ X 



where F s (t) and Fjj(t) are modified noise forces with 
correlations 



V 7 
7i 



(Fl{t)F D {t')) = 2 1 ^ 

(F Sl (t)H 1 ( t ')) = 2 ^- 

1 

(Fl(t)F Sl (t'))=2 lL ^ 



(D12) 

(D13) 

(D14) 

(D15) 
(D16) 



and all other correlations can be neglected. The coher- 
ent part of the interaction can thus be obtained from an 
effective hamiltonian 



H, 



eff 



SíPd + h.c. 



(D17) 



where 



the 



interaction 



rate 



is 



/(NA). 



We note ( [Dll ) that cavity losses are strongly sup- 
pressed in the limit r\ ^> 1. Indeed, subsequent to the 
large group velocity reduction p5[ , the polariton is al- 
most purely atòmic and the excitation leaks very slowly 
out of the medium. The equation of motion for coherence 
(DIC) contains a loss term (due to isotropic sponta- 



neous emission) and a linear gain term (due to emission 
into bright polariton). The two can compensate each 
other. However the linear phase-insensitive amplification 
is also accompanied by correspondingly increased fluctua- 
tions, represented by new Langevin forces Fo(t), Fg (t). 
In the case that g\ — g2 and when all Rabi freqüèn- 
cies are taken to be real, we have the interaction rate 
e = x/V^ = "i^2/A. 
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